AD-A038  625 


UNCLASSIFIED 


GENERAL  APPLIED  SCIENCE  LABS  INC  WESTBURY  N Y F/6  1/3 

THE  EFFECTS  OF  VISCOUS  INTERACTION  ON  THE  TRANSONIC  JET-FLAP. (U) 

HAR  77  P BARONTI*  S ELZWEIG*  G MILLER  F44620-76-C-0017 


GASL-TR-237 


AFOSR-TR-77-0446 


I 


t- 1 — 


\ / 


r-. 

'? 

' ; i 


« : 


/ 


£ 

s 


in 


n i • ' • ' 


General  Applied  Science  Laboratories  inc. 


Approved  for  public  releti3i| 

distribution  unlimited. 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  (AFSC) 
NOTICE  OF  TRANSMITTAL  TO  DDC 
This  technical  report  has  been  reviewed  and  is 
■improved  Tor  publ  e release  I AW  APR  190-12  (7b). 
I-3t.ribution  is  unlimited. 

A.  D.  BLOSE 

Technical  Information.  Qffiaar 


I 

1 

1 

] 

] 

3 


a 


. 


i 

* 


i 


MARCH  1977 


GASL  TR  237 

THE  EFFECTS  OF  VISCOUS  INTERACTION  ON 
THE  TRANSONIC  JET-FLAP 

BY 

P.  Baronti,  S.  Elzweig  & G.  Miller 


A FINAL  REPORT 
PREPARED  FOR 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
BOLLING  AIR  FORCE  BASE 
WASHINGTON,  D.  C.  20332 


UNDER 

CONTRACT  NO.  FAA620-76-C-001 7 


r/VT.VSU'H  N STATEMENT  A 
...ibiic  rolea*ej 
Distri'  Unlimited 


GENERAL  APPLIED  SCIENCE  LABORATORIES,  INC. 

Merrick  and  Stewart  Avenues 
Westbury,  New  York  11590 

1 


'V  a 


I'TTfT'W 


TR  237 
ABSTRACT 


The  inviscid  and  viscous  effects  associated  with  a jet-flapped  airfoil  in 
the  transonic  regime  are  investigated.  The  particular  formulation  developed 
employs  integral  methods  for  the  description  of  the  viscous  portions  of  the 
flow  (utilizing  the  turbulent  kinetic  energy  equation)  and  a relaxation  method 
for  the  inviscid  portions.  The  v i scous- i nvi scid  coupling  procedures  along  the 
wing  and  the  jet  acknowledge  the  importance  and  magnitude  of  the  viscous  ef- 
fects on  pressure  distribution  and  directly  incorporate  the  displacement  ef- 
fect on  the  wing  and  the  effects  of  the  jet  mixing  into  the  boundary  conditions 
for  the  computation  of  the  entire  flow.  Some  numerical  results  are  presented 
to  validate  the  solution  techniques  that  have  been  selected  and  to  assess  the 
importance  of  mass  and  momentum  entrainment  into  the  jet  on  jet  flap  effective- 
ness. 


\ v' 


I 

I 

I 


TR  237 

INDEX 


I.  INTRODUCTION 

II.  INTEGRAL  METHOD  FOR  VISCOUS  JET  SOLUTION 

III.  COUPLING  PROCEDURE  FOR  COMPLETE  SOLUTION 

IV.  RESULTS  AND  CONCLUSIONS 
REFERENCES 

APPENDIX  A 
APPENDIX  B 


1 

4 

9 

20 


-1 L- 


TR  237 

LIST  OF  FIGURES 


FIG.  t. 

THE  VISCOUS  FLOW  FIELD  ALONG  A CURVED  JET  FLAP 

5 

FIG.  2. 

SCHEMATIC  REPRESENTATION  OF  THE  INVISCID  JET-FLAPPED 
AIRFOIL 

10 

FIG.  3. 

SCHEMATIC  OF  THE  E,n  COMPUTATIONAL  PLANE 

13 

FIG.  4. 

EFFECTS  OF  V 1 SCOUS- 1 NV  1 SC  1 D INTERACTION  ON  6% 
CIRCULAR-ARC  AIRFOIL 

21 

FIG.  5. 

EFFECT  OF  JET  MIXING  ON  PRESSURE  DISTRIBUTION  OVER 
A 6%  CIRCULAR-ARC  AIRFOIL 

22 

FIG.  6. 

EFFECT  OF  JET  VISCOSITY  ON  PRESSURE  DISTRIBUTION 
A 12%  CIRCULAR-ARC  AIRFOIL 

OVER 

23 

FIG.  7. 

PRESSURES  AND  NORMAL  VELOCITIES  ALONG  BOUNDARIES 
JET  MIXING  REGION 

OF 

24 

i 


I 

I 

I 

I 


1 


a 

a. ,G,L 
c 


q 

R 

s,n 


u,v 


U 


oo 


x,y 

a 

6 

r 

6 

6* 

e 

V 


S.n 


TR  237 

LIST  OF  SYMBOLS 

speed  of  sound 

Bradshaw  turbulent  parameters 

airfoil  chord  . . 

skin  friction  coefficient 

pressure  coefficient 

jet  momentum  coefficient 

airfoil  shape 

free  stream  Mach  number 

static  pressure 

wake  function 

total  velocity 

local  radius  of  curvature  of  jet  flap 
natural  coordinates 
velocity  components 
free  stream  velocity 
Cartesian  coordinates 

angle  of  attack  of  wing,  also  stretching  parameter 
functional  form,  Equation  (19) 
ci rculat ion 

boundary  layer  thickness 

boundary  layer  displacement  thickness 

angle  related  to  jet  curvature,  Figure  (2) 

relaxation  step 

stretched  coordinates 


it 


Coles  wake  parameter 


I 

I 


TR  237 

LIST  OF  SYMBOLS  (Continued) 


dens i ty 

initial  jet  flap  angle,  also  shear  stress 
perturbation  potential 


stream  function 


SUBSCRIPTS  AND  SUPERSCRIPTS 


e condition  at  boundary  layer  edge 

j jet  flap 

l airfoil  lower  surface 

u airfoil  upper  surface 

B lower  side  of  jet  sheet 

T upper  side  of  jet  sheet 

LE  leading  edge 

TE  trailing  edge 

o condition  along  base  line  of  jet 

1,2  conditions  at  upper  and  lower  boundary  of  jet 

transformed  variables 


-v- 


I 


TR  237 
SECTION  I 
INTRODUCTION 


The  use  of  jet  flaps  to  enhance  maneuverability  in  transonic  flight  has  been 
under  scrutiny  in  recent  years.  Wind  tunnel  experiments  on  jet-flapped  air- 
foil s 1,2  ’ 3 ’ 5 have  indicated  that  significant  lift  augmentation  and  improve- 

ment of  the  drag  polar  at  large  lift  can  be  obtained  in  the  transonic  regime; 
at  the  same  time,  retardation  of  boundary  layer  separation  and  of  buffet  onset 
can  be  achieved.  A body  of  theoretical  work  is  also  available  which  provides 
a mathematical  and  numerical  framework  for  the  solution  of  the  transonic  inviscid 
flow  field  around  jet-flapped  airfoils.  Ives  and  Melnick6  have  numerically 
solved  the  full  transonic  potential  equation  while  allowing  the  jet  flap  to 
maintain  a constant,  finite,  thickness.  Malmuth  and  Murphy  have  given  a 
linearized  solution  for  two-dimensional  and  three-dimensional  jet-flapped  wings 
with  infinitely  thin  jets.  However,  none  of  these  analytical  models  has  ex- 
plored the  viscous  effects  of  the  jet  on  airfoil  performance,  since  jet  mixing 
has  been  excluded  and  the  jet  momentum  assumed  constant  along  the  jet.  As 
originally  indicated  by  Stratford,9  the  effect  of  entrainment  is  to  cause  the 
jet  to  act  like  a sink  with  respect  to  the  inviscid  field,  with  consequent  loss 
of  recoverable  jet  thrust.  More  recently,  Yoshihara  and  Zonars10  have  reviewed, 
although  qual ! tatively,  the  effects  of  jet  viscosity  on  jet  flap  performance  in 
the  transonic  regime. - They-  have -pointed  out,  on  the  basis  of  experimental  evi- 
dence, several  important  phenomena:  the  ejector  action  of  the  jet  in  increasing 

trailing  edge  suction  (and  thus  decreasing  thrust  recovery),  the  effective 
reduction  of  the  jet  exit  angle  because  of  jet  mixing^and  the  importance  of  jet 
and  wake  vi scous- invi sci d interaction  in  providing  a realistic  picture  of  the 
transonic  flow  field  for  large  jet  deflection  angles. 


For  these  reasons,  an  effort  to  assess  the  effects  of  vi scous- invi sci d inter- 
action on  jet  effectiveness  was  initated  by  the  authors  in  Reference  (11).  An 
invisci^  analysis,  as  well  as  an  attendant  computer  program  similar  to  that  of 
Ives  and  Melnick,  was  developed.  In  addition,  integral  techniques  for  the 
separate  analyses  of  the  viscous  regions  along  the  airfoil  and  the  jet  were 
developed  and  the  coupling  procedures  outlined.  The  choice  of  an  integral 
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method,  rather  than  a finite  difference  technique,  was  prompted  by  previous 
work  which  showed  that  boundary  layers  over  supercritical  wings  and  wakes 
can  be  well  treated  by  integral  methods,  once  implemented  by  a solution  of  the 
turbulent  kinetic  energy  equation.  Also,  the  separation  of  the  flow  into 
inviscid  and  viscous  regions  greatly  reduces  the  numerical/computational  effort 
as  compared  to  time-dependent  solutions  of  the  "turbulent"  Navier-Stokes 
equations,  which  may  be  foreseeably  attempted  in  the  future. 

In  the  present  report,  the  integral  solution  of  the  viscous  flow  along  the  jet 
and  the  viscous- invi sci d coupling  along  the  jet  have  been  reconsidered.  Further 
numerical  experimentation  has  in  fact  revealed  that  the  integral  solution  of  the 
viscous  flow  along  the  jet  developed  in  Reference  (11)  was  unstable  when  the 
jet  was  highly  curved.  This  has  necessitated  the  reformulation  of  the  momentum 
integral  equations  in  terms  of  velocity  profiles  expressed  by  polynomials  of 
second  degree  (instead  of  third  degree,  as  it  was  originally  proposed  in  Refer- 
ence (11)).  In  addition,  a more  consistent  vi scous- invi sc i d coupling  procedure 
along  the  jet  has  been  elaborated. 

The  present  report  completes  the  formulation  of  the  problem  of  viscous- i nvi sci d 
interaction  along  jet-flapped  airfoils,  initiated  in  Reference  (11).  The 
salient  features  of  the  method  are  reviewed  here.  They  consist  of  an  inviscid 
relaxation  procedure  that  resembles  that  of  Ives  and  Melnick,  integral  methods 
for  the  description  of  the  viscous  (turbul ent) portions  of  the  flow  along  the 
airfoil  and  along  the  curved  jet,  and  an  appropriate  coupling  procedure  for  the 
treatment  of  the  vi scous- i nvi sc i d interaction  along  the  wing  and  along  the  jet. 

The  integral  solution  of  the  jet  flap  presented  here  suffers  from  some  limita- 
tions, since  it  requires  a modelling  of  the  viscous  jet  flap  which  may  not  always 
be  consistent  with  physical  reality.  In  fact,  the  present  representation  of  the 
velocity  field  at  the  wing  trailing  edge  or  equivalently  at  the  jet  exit,  con- 
siders only  the  jet  velocity  profile  and  disregards  the  detailed  velocity  pro- 
files of  the  boundary  layers  on  both  sides  of  the  wing  which  merge  with  the 
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jet  flow  (however,  in  the  actual  computations  performed  here,  the  initial  jet 
profile  is  modified  so  as  to  satisfy,  on  the  average,  the  mass  and  momentum 
flux  of  the  jet-boundary  layer  combination).  This  representation  is  realistic 
when  the  boundary  layer  on  the  wing  is  thin  and  attached,  or  when  it  has  very 
small  separated  regions  at  the  wing  trailing  edge,  but  it  becomes  invalid  in 
the  case  of  jets  with  large  initi.al  deflection  angles  which  may  induce  boundary 
layer  separation  and  large  recirculation  regions  at  the  wing  trailing  edge.'^ 
Notwithstanding  these  limitations  it  appears,  however,  that  the  integral  repre- 
sentation may  be  usefully  and  simply  invoked  for  a quantitative  determination 
of  viscous  effects  on  jet-flap  performance  in  many  configurations  of  practical 
interest. 


The  report  is  organized  as  follows.  In  Section  II  the  integral  analysis  of  the 
viscous  region  along  a curved  jet  is  presented.  Section  111  outlines  the 
inviscid  analysis  and  the  coupling  procedure  of  the  viscous  and  inviscid  solu- 
tions along  the  wing  and  along  the  jet.  Some  numerical  results,  indicating 
the  effects  of  jet-mixing  and  entrainment  on  the  pressure  distribution  over  a 
two-dimensional  airfoil,  as  well  as  some  conclusions  with  respect  to  the  quanti- 
tative effect  of  jet  momentum  and  deflection  angle,  are  presented  in  Section  IV. 
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SECTION  I I 

INTEGRAL  METHOD  FOR  VISCOUS  JET  SOLUTION 


The  authors  have  developed  in  Reference  ( 1 1 ) an  integral  solution  of  the  vis- 
cous flow  along  a curved  jet.  An  integral  solution  of  the  boundary  layer  flow 
over  an  airfoil,  under  subcritical  and  supercritical  conditions,  had  been 
previously  developed  in  Reference  (12).  Both  solutions  utilize  an  integral  so- 
lution of  the  turbulent  kinetic  energy  equation13  and  simple  coordinate  trans- 
formations to  account  for  compressibility  effects.  The  formulation  of  the  in- 
tegral solution  of  the  viscous  region  along  a curved  jet  is  now  reconsidered. 


With  reference  to  Figure  (l),  a jet  is  issued  from  the  trailing  edge  of  an  air- 
foil with  a deflection  angle  t.  The  flow  around  the  wing  is  assumed  to  be  in- 
viscid  (thus  the  boundary  layer  over  the  wing  is  neglected)  and  only  the  mixing 
of  the  jet  with  the  inviscid  flow  field  is  considered.  A line  s,  not  neces- 
sarily the  centerline  of  the  mixing  region,  is  selected  as  a "base"  streamline 
of  the  flow.  Its  position,  to  be  determined  by  the  viscous-invi scid  coupling 
discussed  below,  can  be  proved14  to  be  accurate  to  order  6 of  boundary  layer 
theory.  Along  this  streamline  v = 0,  but,  still  within  the  accuracy  of  bound- 
ary layer  theory,  the  velocity  profile  u(n)  may  have  a slope  different  from 
zero.  Thus,  a shear  different  from  zero  is  allowed  along  the  s-st reaml i ne . 

The  mixing  boundaries  of  the  curved  jet,  defined  by  the  boundary  lines  n=S. 
and  n=$2  , are  asymmetric  with  respect  to  the  "base"  streamline  and  the  invis- 
cid velocities  u^  and  u^  along  6^  and  6^  are,  in  general,  different  from  each 
other.  The  problem  solution  requires  a representation  of  velocity  and  shear 
profiles  which  accounts  for  these  asymmetries. 


In  Reference  (11),  the  velocity  profiles  were  represented  by  the  third  degree 
po 1 ynom i a 1 


u = u„  + a.  n + a. 


+ a. 


with  the  coefficients  a^  and  a^  defined  by  satisfying  the  boundary  conditions 
u = Uj  at  n = Sj  and  u = u^  at  n = 6^  (note  that  the  edge  derivatives  are  left 
unspecified).  The  flow  development  was  then  defined  in  terms  of  the  dependent 
variable  Ug(s),  6^  (s)  , 62(5)  and  a^(s).  These  were  obtained  by  integrating  the 
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streamwise  momentum  equation  four  times,  that  is  from  the  "base"  streamline 
n = 0 to  6^/2  and  to  6^  and  from  n = 0 to  6^/2  and  to  In  a similar  fashion, 

in  order  to  obtain  the  shear  distribution  across  the  boundary  layer  (needed  in 
the  integration  of  the  momentum  equation),  shear  profiles  were  prescribed  by  the 
fourth  order  polynomial 

2 3 A 

t = b^  + b1  n + b^  n + b^  n + b^  n 

with  the  coefficients  b^ , b^  and  b^  obtained  by  satisfying  the  conditions  x = 0 
at  n = and  n = 6^  and  by  satisfying  the  s-momentum  equation  along  the  base 
streamline.  The  coefficients  b^  and  b^  were  then  obtained  by  integrating  twice 
the  turbulent  kinetic  energy  equation  of  Bradshaw,  Ferris  and  Atwell be- 
tween n = 0 and  6^  and  between  n = 0 and  6^- 

The  resulting  system  of  integral  momentum  equations  was  analyzed  extensively  for 
asymmetric  jet  flows  but  was  found  to  be  unstable  for  cases  of  large  deflection 
angles  and  large  differences  between  u^  (s)  and  U2(s).  The  instabilities  were 
traced  to  the  representation  of  the  velocity  profiles  by  polynomials  of  the 
third  degree,  which  could  produce  negative  velocities  and  uncontrolled  velocity 
ampl if icat ions.  A representat ion  of  velocity  and  shear  profiles  by  lower  degree 
polynomials  has  then  been  investigated.  Resort  has  been  made  to  parabolic  velo- 
city profiles  and  third  degree  polynomials  for  the  shear.  An  approach  which 
utilizes  the  integration  of  the  momentum  equation  in  the  three  steps  from  0 to 
6^/2  , from  0 to  6^/2  and  from  6^  to  was  then  established  and  coupled  to  an 
integral  solution  of  the  turbulent  kinetic  energy  equation  from  6^  to  62-  This 
system  was  found  to  be  stable  even  for  the  case  of  severe  asymmetries  between 
the  velocities  u^  (s)  and  U2(s).  It  was  also  found  to  be  as  accurate  as  the 
previous  system  based  on  polynomials  of  third  degree  for  the  velocity,  in  those 
cases  where  the  latter  system  gave  stable  solutions. 

The  details  of  the  formulation  are  as  follows.  A restricted  Dorodn i tsyn-Howarth 
compressibility  transformation  (where  only  the  n coordinate  is  stretched)  is 
first  invoked,  see  Appendix  A,  to  yield  the  transformed  momentum  equation  (bar- 
red quantities  identify  the  transformed  variables) 
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(1) 


9u  3u 

u — — + v — — 

9 s 9n 


. 1 i£  + l il 
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The  velocity  profiles  are  described  by  a second  degree  polynomial 


-2 


u = u + a.n  + a„n 
o 1 2 


and  the  shear  profile  by  a third  degree  polynomial 


(2) 


T = bg  + b^n  + b^ (3) 

The  coefficients  a1  and  a^  in  Equation  (2)  are  defined  by  satisfying  the  bound- 
ary conditions  u = u^  at  n = 6^  and  u = u^  at  n = 8^.  The  coefficients  b^  , b^ 
and  bg  of  Equation  (3)  are  obtained  by  satisfying  the  conditions  t = 0 at 
n = and  n = 6^  and  by  satisfying  the  s-momentum  equation  along  the  base 
streamline. 


A Crocco  integral,  see  Appendix  A,  for  the  definition  of  the  density  varia- 
tion across  the  viscous  layer  is  then  assumed  to  be  locally  valid.  In  addi- 
tion, a linear  variation  of  the  pressure  across  the  jet  can  be  proved,  see 
Appendix  B,  to  be  adequate  even  in  the  case  of  highly  curved  jets.  The  inte- 
gration of  the  momentum  equation.  Equation  (l),  between  0 and  6^/2,  0 and 
^2/2  and  between  6^  and  6 ^ can  then  be  carried  out,  see  Appendix  B,  to  yield 
the  system  of  equations 


A11  A12  A1 3 

duQ/ds 

C1 

A21  A22  A23 

d6^ /ds 

= 

C2 

A,,  A„„  A,, 

d^/ds 

C, 

31  32  33 

3 

in  the  variables  u , 6,  and  6„.  The  coefficients  A.  . are  functions  of 

°_  < 2 i,j 

u0,  6^,  62*  u^  and  U2,  The  coefficents  C.  are  functions  also  of  the  values  of 
the  shear  at  n = 0,  6^/2  and  62/2,  and  thus  of  the  coefficient  b^  of  Equation 
(3).  The  coefficient  b^  is  obtained  by  integrating  the  turbulent  kinetic 
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energy,  TKE,  equation  across  the  jet. 

It  was  shown  in  Reference  (12)  that  the  compressible  TKE  equation  of  Bradshaw 
and  Ferris16  can  be  transformed  by  the  Dorodn i tsyn-Howarth  transformation  into 
the  incompressible  form  given  in  Reference  (15)  and  that  such  an  equation,  de- 
rived originally  for  boundary  layer  flows,  can  be  also  applied  to  wakes.  The 
applicability  of  this  equation  to  straight  and  curved  jet  is  now  assumed.  Thus, 
integration  between  <5^  and  6^  °f  the  equation 


- 1/2  _ . . 3/2 

u 4 (-4-)  + v 4 (-4)  - 1 4 + (44)  4 (G  D + 44—  = 0 (5) 

3s  2a. p 3n  2a. p p 3n  p 3n  p C 


3n  2a. p 
1 


yields,  as  shown  in  Appendix  B,  an  equation  for  b^  of  the  form 


db 

Bi  *r=  Gi 


Equations  (4)  and  (6)  provide  the  closure  of  the  system  of  equations  and  un- 
knowns. Application  of  the  transformation  in  an  inverse  manner  provides, 
finally,  the  desired  flow  quantities  in  the  physical  plane. 
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SECTION  I I I 

COUPLING  PROCEDURE  FOR  COMPLETE  SOLUTION 


A vi scous- i nvi sc  id  coupling  procedure  which  accounts  for  strong  interaction 
along  the  wing  and  the  jet  is  now  developed.  The  boundary  layer  solution 
along  the  wing  is  the  one  developed  in  Reference  ( 1 2 ) . The  viscous  solution 
along  the  jet  was  presented  in  the  preceding  section.  The  details  of 
the  initial  cond i t ions  at  the  wing  trailing  edge,  where  the  boundary  layer 
flow  merges  with  the  jet  flow,  are  disregarded.  The  occurance  of  separation 
and  of  recirculation  regions  at  the  trailing  edge  are  excluded;  however,  ini- 
tial jet  profiles  are  chosen  such  that  the  mass  and  momentum  flux  of  the  jet- 
boundary layer  combination  is  conserved.  The  viscous  analyses  are  coupled 
with  a relaxation  sol ut ion  of  the  external  inviscid  flow.  The  inviscid  solu- 
tion used  here  resembles  that  of  Ives  and  Melnick6  and  is  briefly  reviewed  in 
order  to  introduce,  in  a sequential  manner,  the  description  of  the  coupling. 

A schematic  representation  of  the  inviscid  airfoil/jet  flap  combination  with 
related  nomenclature  is  shown  in  Figure  (2).  Solutions  for  the  flow  field  in 
the  inviscid  limit  are  sought  by  numerical  integration  of  the  full  isentropic, 
compressible  flow  equations  written  in  terms  of  the  perturbation  potential 
4> * scaled  by  the  freestream  velocity  and  chord,  and  arranged  in  the  canonical 
form 

(1  - q2/a2)  $ss  + $nn  = 0 (7) 

where  s and  n are  othogonal  streamline  coordinates,  a is  the  local  sound 
speed  and  q is  the  total  velocity. 


For  the  inviscid  problem  the  boundary  conditions  along  the  airfoil  and  along 
the  jet  are,  respectively,  (it  will  be  seen  below  that  the  coupling  procedure 
will  manifest  itself  only  in  a change  in  boundary  conditions) 


V(1  4 *x)]*,u  = 


(8) 
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and 


+ 4>„)  1 = Y;(x) 

y A j j 


(9) 


AJso,  the  pressure  difference  across  the  jet  is  required  to  satisfy  the 
centrifugal  force  balance 


Ac  = - c./R 
P J 


(10) 


where  c.  is  the  jet  momentum  per  unit  width  (divided  by  the  free  stream 

J2 

OO  CIO 


dynamic  pressure  1/2  p^U  times  the  airfoil  chord  c)  and  R is  the  dimension- 
less (scaled  to  c)  radius  of  curvature  of  the  jet.  Finally,  boundary  condi- 
tions are  imposed  at  infinity.  For  < 1 , when  the  problem  is  inherently 
elliptic,  the  standard  linearized  vortex  solution  is  imposed  on  the  far  field 
potential , viz. 


<ra+ri)  , -1 

- — tan 


/rrj 


2it 


(11) 


where  T is  the  circulation  about  the  airfoil  and  T.  is  the  circulation  a- 
a J 


bout  the  jet  sheet. 


To  numerically  solve  the  system  (7)  - (11),  a transformation  from  physical 
space  into  a unit  square  is  performed.  The  mapping  functions  are 


€ = tanh  (cijx) 


(12) 


n = 


where 


g(x) 


(a2 

[y  - 

g(x)] } 

yL.E. 

for 

x < xL  E 

fu(x) 

for 

XL.E.  < X < XT . E . 1 

r y > fu(x) 

fA(x) 

for 

XL.E.  < X < XT.E. 

. y < ^(x) 

Yj(x) 

for 

X " XT.E. 
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A schematic  representation  of  the  computat ional  plane  is  shown  in  Figure  (3). 

The  solution  of  Equation  (7)  is  performed  in  £,n  using  the  rotated  difference 
technique  of  Jameson17.  Tridiagonal  equations  of  the  form 


a<J>.  . . + b<f> . + cij>.  ...  = d 


i ,j~1 


. + c <(>.  . 

,J  i»J  + 


(13) 


are  derived,  where  a,  b,  c and  d are  functions  of  the  coordinates  and  the 
previous  iteration  values. 


The  jet  boundary  condition,  Equation  (10),  is  incorporated  into  the  difference 
scheme  following  the  approach  of  Ives  and  Melnick.  Across  the  jet 


4>g  + r(s) 

where  the  subscripts  T and  B denote  the  top  and  bottom  of  the  jet  sheet,  re- 
spectively, and  T denotes  the  circulation  at  a given  streamwise  station  along 
the  jet.  Since 


(1M 


dr 

ds 


:nr  = 
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I 

"i 

*1 

-fi  i 

l!  1 


i 

l 
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t I 


and 


ds  = Rde 


Equation  (14)  can  be  combined  with  Equation  (10)  and  integrated  to  yield 


+ C 


0 


qT  “ qB 


PB  " PT 


TE 


de 


(15) 


Accordingly,  the  finite  difference  approximation  to  Equation  (7)  just  below 
and  just  above  the  jet  becomes 
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"B  ♦i.J-l  + bB  ♦l,jB  + CB  ♦ilj+1B  " dB 


’t 


+ bT  *i,j  + Ct  ♦:  : j-i  = d 


T Ti,j+1  T 


(16) 


where  . . and  4>  . . are  fictitious  points  outside  the  flow  field,  and 

1 » J 1 g ■ > J ~ 1 y 

where  <J>.  . and  4> . . are  the  values  of  the  potential  just  below  and  just  a- 
1 * J B 1 » JT 

bove  the  jet,  respectively. 


The  requirement  that  the  velocity  normal  to  the  jet  sheet  is  continuous  across 
the  sheet  yields  an  additional  relation  which,  once  combined  with  Equations 
05)  and  (16),  eliminates  the  fictitious  potentials  <f>.  . and  Thus, 


the  difference  formulation 


J + 1 


B 


i , j -1. 


b„  a. 


3T  0 + ^7}  *i,j-1  + (bT  +~^-]  ♦i.J  + (aT  +CT}  *|, 


j + 1 


(17) 


aT_d 
'B 


I— ^ + dT  - bT  r + aT  ^ ^ 0 - S2)  Sin26  (^  - «U 

2 q. 


valid  along  the  jet  flap  is  achieved.  The  jet  shape  is  obtained  by  tracing 
the  streamline  coming  off  the  airfoil  trailing  edge  and  is  updated  during 
each  interaction. 


The  coupling  of  the  viscous  results  to  the  inviscid  analysis  is  accomplished 
by  a modification  of  the  boundary  conditions  along  the  airfoil  (Equation  8)  and 
the  jet  (Equations  10  and  17).  For  the  flow  over  the  wing,  a strong  interac- 
tion coupling  is  utilized.  The  integral  of  the  continuity  equation  in  the  vis- 
cous region  yields 


d6 

dx 


. d tn(p  u ) 
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In  addition,  as  shown  in  Reference  (11),  one  can  express  the  displacement 
distribution  in  the  form 


+ 


d ln(p  u ) 
e e 

dx 


(19) 


(where  BQ  and  Bj  are  functions  of  the  boundary  layer  quantities  c^,  ir,  6 
for  attached  flow  and  P and  S*  for  separated  flows,  see  References  11  and  12). 
Combining  Equations  (18)  and  (19)  yields 


v 

e 

u 

e 


d 2.n(p  u ) 
e e 

dx 


where  B2  = 8q 


(6  - 6 ) . Thus,  the  inviscid  flow  boundary  condition 


$y  - f (x) 


on  the  surface  of  the  airfoil  with  profile  y = f(x)  is  replaced  by 


<t>y  = f (x)  + 61  + B3  4>xX,  (2°) 

2 * 

where  B,  = 8~  (1  + M ),  along  the  line  y = f(x)  +6  (x) , or  by  the  condi- 
i 2.  <*> 

t ion 


+ 


8 


v-1 

1 


where  the  superscript  v denotes  the  current  iteration  step.  The  coupling  of 
inviscid  and  viscous  flows  is  provided  by  the  coefficients  B^  and  and  by 
the  distribution  6 (x)  (through  the  positioning  of  the  boundary  condition) 
which  are  functions  of  the  boundary  layer  properties  (but  not  of  their  grad- 
ients). This  method  of  coupling  is  similar  in  philosphy  to  that  of  Tai18  in 

that  the  viscous  effect  on  displacement  (through  the  variable  v ) is  directly 

e 

incorporated  into  the  solution  for  the  inviscid  field.  It  is  also  noted  that 
in  the  "full"  solution  of  the  Nav ier-Stokes  equations,  Diewart19  utilizes  the 
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flow  velocity  v^  at  the  edge  of  the  "viscous"  layer  directly  to  determine 
the  outer  flow  properties. 

With  regard  to  the  vi scous- inv i sc  id  interaction  along  the  jet,  the  two  effects 
of  momentum  and  mass  entrainment  must  be  recognized  and  incorporated  directly 
into  the  coupling  procedure.  The  first  modifies  the  definition  of  c..  The 
second  is  the  sink  effect  and  implies  a discontinuity  of  normal  velocity  a- 
cross  the  viscous  jet.  To  determine  the  effect  of  entrainment  on  the  jet 
momentum  consider  the  equation  for  pressure  across  the  jet 

|£-  = - pu2/R 


and  integrate  it  (in  inner  viscous  variables)  to  yield 


Ap 


pu 


dn 


But,  from  the  outer  (inviscid)  solution 


P = p,  + (t^~)  n + — 
“>  1 3n  j 00 


P-»=  P2  + (^2  + — 


where 


■ • Pr?/R 


/ 3p\  _ 2 . 

(3^2  ~ " P2U2/R 


Thus,  the  effect  of  mixing  on  the  pressure  difference  across  the  jet  can  be 
expressed  by  the  relationship 
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Ac  = - - 
P R 


(pu2  - p2u2) 

p jin 


dn  + 


(pu2  - P^2) 
P U2/2 


dn 


2 2 2 

which  replaces  Equation  (10).  Note  that  for  pu  >>  P^u^>  P2U2’  e*Pres" 

sion  reduces  in  the  inviscid  limit  to  Ac  =-c./R.  Then,  if  an  "effective" 

- . . P J 
momentum  coefficient  c^  is  defined,  viz. 


r(s,- 


(pu2  - P2u2^ 
P U2/2 

OO  CO 


dn  + 


J. 


(pu2  - p^2) 
p U2/2 


dn, 


Equation  (15)  can  be  replaced  by 


J 


= r 


T.E. 


T.E. 


v-1 


^7  (5)  (PT  - qB} 


(c 


v ^ 
c ) 
Pt 


de 


(21) 


to  account  for  momentum  entrainment. 

The  sink  effect  is  accounted  for  by  imposing  along  the  jet  baseline  the 
normal  velocity  discontinuity 

(Av)  = v^"1  - v"'1  (22) 

62  1 

where  v<  and  v,  are  obtained  by  integrating  the  continuity  equation  from 
°2  °1 

the  base  streamline.  Thus 


where  I\  and  Av  are  given  by  Equations  ( 2 1 ) and  (22),  respectively,  in  terns  of 
va lues  establ i shed  at  a preceding  viscous  iteration  step. 


The  replacement  of  Equations  (8)  and  (17)  with  Equations  (20)  and  (2b)  con- 
stitutes the  entire  change  in  the  inviscid  program  to  account  for  strong 
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viscous  interaction.  The  effects  of  viscosity  are  contained  in  the  coef- 
ficient Bj  and  of  Equation  (20)  and  the  coefficients  c^  and  Av  of  Equa- 
tions (21)  and  (22).  The  values  of  these  coefficients  are  obtained  from  the 
viscous  analysis  along  the  airfoil  and  along  the  jet,  developed  in  Reference 
(12)  and  in  the  present  report. 
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SECTION  IV 

RESULTS  AND  CONCLUSIONS 


In  Reference  ( 1 1 ) , preliminary  results  were  reported  to  ascertain  the  importance 
of  viscous  effects  in  the  analysis  of  jet-flapped  airfoils  in  the  transonic  re- 
gime. For  the  airfoil  without  a jet  flap  the  viscous  interaction  model  along 
the  airfoil  was  validated  by  the  results  shown  in  Figure  (b ) relative  to  a 6% 
circular-arc  airfoil  at  Mm  = 0.908.  The  numerical  result  agree  quite  favor- 
ably with  the  experiment  of  Reference  (20) . The  influence  of  jet  mixing  on  the 
pressure  distribution  over  the  wing  can  be  of  equal  importance,  as  manifested 
by  the  result  of  Figure  (5)  which  presents  the  pressure  distribution  over  an 
airfoil  at  zero  angle  of  attack  without  and  with  a straight  jet  flap.  Additional 
results  have  now  been  obtained  by  use  of  the  analysis  developed  in  the  previous 
sections. 


The  influence  of  mass  and  momentum  entrained  by  the  jet  on  the  pressure  distri- 
bution over  an  airfoil  was  studied  by  examining  the  flow  about  a 12%  thick 
circular-arc  airfoil  (which  is  similar  to  a NACA  0012  except  at  the  nose)  at 
zero  angle  of  attack  with  a jet  flap  having  an  initial  deflection  angle  of  30°, 

a jet-momentum  coefficient  c.  = 0.1  and  a mass  flow  rate  p.u.S./p  U c = 0.0b, 

J J J J " “ 

at  a free  stream  Mach  number  of  0.7.  Figure  (6)  presents  the  pressure  distri- 
bution over  the  airfoil  for  the  case  of  an  inviscid,  infinitely  thin  jet  flap 
and  of  a viscous  jet  flap.  For  the  viscous  jet,  the  pressure  distribution  along 
the  upper  and  lower  jet  boundaries  are  given  in  Figure  (7).  Also  shown  in  the 
figure  are  the  corresponding  components  of  the  normal  velocity  as  computed  by 
Equation  (23)  of  the  viscous  analysis.  Note  that  for  this  case  (a  moderate  c^ 
and  a moderate  deflection  angle),  the  sink  effect  near  the  trailing  edge 
changes  the  initial  deflection  at  ole  by  about  6%  (since  tan  30°  = 0.577  and 

(— ) 0.03 b) . Note  also  that  this  effects  damps  quickly  to  extremely  small 

u max 

values  in  a distance  of  about  20%  of  the  chord. 


The  results  of  Figure  (6)  indicate  that  the  pressure  distribution  on  the  top- 
aft  section  of  the  airfoil  is  slightly  altered,  while  that  on  the  bottom-aft 
section  is  essentially  unaffected.  The  lower  pressure  on  the  aft-top  surface 
results  in  a decrease  in  thrust,  as  compared  with  the  inviscid  case,  while  the 
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lift  is  somewhat  enhanced.  The  calculation  also  indicates  that,  for  this  parti- 
cular case,  the  effect  of  momentum  entrainment  by  the  jet  is  small  and,  to  first 
order,  the  jet  momentum  as  a function  of  position  (e.g.  the  Cj(s)  coefficient 
*n  Equation  (2 l) ) i s approximatley  equal  to  the  initial  jet  momentum  c^  . This 
fact  may  not  be  true  for  much  smaller  values  of  Cy  i.e.  for  Cj  <<  0.1.  The 
major  effect  of  entrainment  is  due  to  the  sink  effect  (the  parameter  Av  of 
Equat ion  (22) ) whose  magnitude  can  be  ascertained  from  Figure  (7).  The  effect  of 
jet  entrainment  is  more  pronounced  on  the  airfoil  upper  surface  than  on  the  low- 
er surface,  since  the  flow  around  the  lower  surface  must  experience  a large  pres- 
sure increase  near  the  trailing  edge  because  of  the  initial  jet  deflection, 
while  the  upper  surface  exhibits  a relief  on  adverse  gradient  due  to  the  jet. 

The  top  surface  is  also  more  susceptible  to  the  ejector  effect  of  the  jet  flap. 
In  addition,  the  sink  effect  is  proportional  to  the  magnitude  of  c^  as  indicated 
by  numerical  experiments  carried  with  the  present  analysis.  Thus,  for  a given 
initial  jet  deflection  angle  the  effect  of  viscosity,  through  mass  entrainment, 
increases  with  increasing  jet  momentum. 


The  other  parameter  whose  variation  can  significantly  alter  the  effect  of  en- 
trainment is  the  initial  jet  deflection  angle  t.  As  t increases,  the  inviscid 
effects  of  rear  stagnation  on  the  bottom  surface  and  of  pressure  relief  on  the 
top  surface  become  more  and  more  pronounced.  At  the  same  time,  the  effect  of 
the  jet  on  separation  and  bubble  development  is  enhanced.  Thus,  as  pointed 
out  in  Reference  ( 1 0) , the  identity  of  the  jet  can  be  obliterated  by  viscous 
mixing.  It  therefore  seems  that  the  effect  of  viscosity  is  important  for  small 
(see  Figure  5)  and  large  (see  Reference  10)  jet  deflection  angles.  Thus,  one 
can  expect  the  effects  of  entrainment  will  be  most  significant  for  large  c^. 
and  for  both  small  and  large  t.  Therefore,  a program  which  computes  the 
vi scous- i nvi sci d interaction  is  necessary  to  properly  ascertain  the  lift  and 
thrust  increments  (or  decrements)  under  a particular  set  of  conditions.  The 
coupling  procedure  described  in  the  previous  section  for  the  airfoil  and  the 
jet,  when  utilized  in  concert,  can  adequately  describe  the  effects  of  viscosity 
and  can  be  readily  incorporated  into  existing  codes  for  an  assessment  of  jet 
flap  performance  for  flows  without  large  separation.  In  addition,  as  pointed 
out  in  Reference  ( 1 0 ) , a proper  assessment  of  viscous  effects  for  the  treatment 
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of  jet-flapped  wings  of  finite  span  would  entail  the  development  of  a three- 
dimensional  analysis,  perhaps  utilizing  some  of  the  concepts  developed  herein. 
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APPENDIX  A 

COMPRESSIBILITY  TRANSFORMATIONS  AND  CROCCO 
RELATION  FOR  VISCOUS  JET  ANALYSIS 


A restricted  Dorodni tsyn-Howarth  transformation 


s = s 


n 

n - j 


p/p  dn 


* = <P 

with  barred  quantities  denoting  transformed  values,  is  utilized  here  to 
account  for  compressibility  effects  along  the  jet. 

The  transformation  implies  that 


vp  , ~ 3n 
v = _+  u _ 

P 

Then,  if  it  is  further  assumed  that  t = r,  the  momentum  equation  along  the 
jet  becomes  (to  order  <$) 


- 3u  - 3n 
u — — + v — =■ 

3s  3n 


1_  J_  3t_ 

P 3s  p 3n 


(A, 2) 


Since  across  the  jet 


application  of  the  transformation  also  implies 


A! 


iE.  = ie.  = iE. 

3s  3s  3s 


The  density  term  1/p  in  the  momentum  equation  can  be  expressed  directly  in 
terms  of  velocity  and  pressure  profiles  and  of  edge  conditions  by  invoking 
the  local  validity  of  the  Crocco  integral.  One  then  has 
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APPENDIX  B 

INTEGRAL  METHOD  FOR  CURVED  JET 


Integration  of  the  transformed  momentum  equation,  Equation  (A, 2) , between  0 
and  6^/2,  0 and  and  between  6^  and  6 yields  (for  convenience,  all  trans- 

formed quantities  are  now  indicated  by  unbarred  symbol  whereas  p indicates 
the  compressible  density) 


u dn  = u 


6./2  ds 

i 


udn  = 


<5./2 

■f  i 


— If-  dn  + (x/p)  - (x/p). 

P 85  6i/2  0 


where  i = 1,2,  and 


2 . d 

u dn  - u,  — 
1 ds 


udn  + / udn 


/- 


-h 

K p 


|£-  dn 
9s 


The  density  p has  been  defined  above  by  Equation  (A,**)  , the  velocity  profiles 
are  expressed  by  the  polynomial 


u = uQ  + a^  n + a^n 


where  u^  is  the  veloc  i ty  along  the  "base"  streamline  and  where  the  coefficients 
a^  and  a ^ axe  obtained  by  satisfying  the  boundary  conditions  u=Uj  at  n=6^  and 
u=u^  at  n=6 ^ to  yield 
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[(uj  - Uq ) 6^  ~ (u2  - uQ)  6^]/D 


= I(u2  - uQ)  6t  - (u  - uQ)  62]/D 


5,  {2  ' S2 


The  pressure  distribution  across  the  layer  can  be  obtained  by  integrating 
the  normal  momentum  equation  (A, 3)  to  yield  a fifth-order  polynomial  for  p. 
This  was  originally  proposed  in  Reference  (11).  It  was  felt,  however,  that 

since  3p/3n  cannot  change  sign  across  the  layer,  it  may  be  possible  to  lower 
the  order  of  the  polynomial.  Extensive  investigation  has  indicated  that  a 
linear  solution  is  adequate.  A pressure  distribution  of  the  form 


[{-'  -y2)  I(P2  s!  - P1  52>  + (P1  - p2> 


has  been,  henceforth,  utilized  here. 

Thus , by  i nsert i ng  Equations  (A, 4) , (B,3)  and  (B,4)  into  Equations  ( B , 1 ) and 
(B,2)  and  performing  the  integrations  and  differentiations  one  obtains  the 
system  of  equations 


^0  + a _2_ 
J11  ds  a12  ds  13  ds 


^1  ^2  r 

a21  ds  a22  ds  a23  ds  " C2 


(B,5) 


du  d<5  d6 

a,.  0 t 1 . 2 - 

31  ds-  a32  ds“  a33  ds~  = C3 


which  defines  the  flow  variables  u^.S^  and  62> 


To  obtain  the  shear  distribution  and,  thus,  the  values  T62/2  anci  T0 
appearing  in  Equations  ( B , 1 ) and  (B,2)  it  is  assumed  that  (unbarred  quanti- 
ties shall  indicate  transformed  values). 


Li 


B2 
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r/p  = bg  + b^n  + b^n^  + b^n^ 


(B  ,6) 


The  coefficient  b^  is  obtained  by  satisfying  the  s-momentum  equation  along 
the  n = 0 line;  hence 

' " 0 0 ds  ?0  Ss  0 


The  coefficients  b and  b„  are  obtained  by  satisfying  the  boundary  conditions 
0 2 

x = 0 at  y = 61  , 62 


Hence, 


b0  ■ i;  ‘-(bl  ei  * b3  S?>  s2  * <b>  *2  + b3  s2>  si} 


where 


b2  = {"(bl  61  + b3  62}  + (b1  61  + b3  61)} 


2 2 
D2  = 62  - 6, 

The  coefficient  b^  is  then  obtained  by i nteg rat i ng  the  TKE  equation  (Equation  B,5) 
across  the  viscous  jet  layer,  viz., 


1 _d_ 

2a.  ds 


— udn 


B3 
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with  the  parameters  a.  and  L given  in  Reference  (15).  By  inserting  Equa- 
tions (B, 3)  and  (B,6)  into  Equation  (B,7),  on  then  obtains  the  ordinary  dif- 
ferential equation  for  of  the  form 


# 


B 


1 


G 


1 


Thus,  Equations  (B,5)  and  (B,7)  provide  the  sought  closure 
unknowns  for  the  solution  of  ^(s),  62(s),  uQ(s)  and  b^(s) 
plane. 


of  equations  and 
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